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Abstract 

We study Andreev reflection (AR) tunneling through a quantum dot (QD) 
connected to a ferromagnet and a superconductor, in which the intradot 
spin-flip interaction is included. By using the nonequibrium-Green-function 
method, the formula of the linear AR conductance is derived at zero temper- 
ature. It is found that competition between the intradot spin-flip scattering 
and the tunneling coupling to the leads dominantes resonant behaviours of 
the AR conductance versus the gate voltage. A weak spin-flip scattering leads 
to a single peak resonance. However, with the spin-flip scattering strength 
increasing, the AR conductance will develop into a double peak resonannce 
implying a novel structure in the tunneling spectrum of the AR conductance. 
Besides, the effect of the spin-dependent tunneling couplings, the matching 
of Fermi velocity, and the spin polarization of the ferromagnet on the AR 
conductance is eximined in detail. 
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1. Introduction 

With the advances of nanofabrication and material growth technologies, it has been pos- 
sible to realize various kinds of hybrid mesoscopic structures 1-4 . Recently, spin-dependent 
electronic transport through these hybrid mesoscopic structures has become one of the major 
focuses of the rapidly developing spintronics 5 for both its fundamental physics and potential 
applications. In particular, the Andreev reflection (AR) in spin-polarized transport through 
ferromagnet-superconductor (F-S) junctions has been examined based on the scattering ma- 
trix formulation 6-11 . It is found that in the case of low bias voltage, AR tunneling at the 
F-S interface is strongly affected by the spin polarization of the ferromagnet side 6 and the 
measuring of the differential AR conductance can successfully determine the spin polariza- 
tion at the Fermi energy for several metals 7 . In addition, further calculations 8,9 showed that 
the AR conductance of F-S junction is also modified by the Fermi velocity mismatch, and 
it may even appear a peak with the varying of spin polarization of the ferromagnet. 

On the other hand, spin-dependent resonant tunneling through a quantum dot (QD), a 
small system characterized by discrete electronic states, coupled a feeromagnet (F) and a 
superconductor (S) leads has been another interesting subject of experimental and theoret- 
ical investigations for the past decade. Zhu 12 et.al. proposed an efficient mechanism for the 
operation of writing spin in such the F-QD-S system with the principle of the Andreev re- 
flection induced spin polarization. They 13 also studied the AR tunneling through a QD with 
two ferromagnets and a superconductor, in which only one spin-degenerate level of the QD 
is considered and the intradot Coulomb interaction is ignored. In this three terminal hybrid 
structure, the transport conducted by crossed AR, which involves an incident electron with 
spin a from one of the ferromagnets picks up another electron with the opposed spin a from 
the other one, both entre the S-lead and form a Cooper pair, is particularly interesting. It is 
found that the AR tunneling processes are, besides the spin polarizations and the matching 
condition of Fermi velocity, strongly dependent on the title angle between the magnetization 
orientations of the two F-leads. Feng and Xiong 14 investigated the transport properties of 
a F-QD-S system, in which both the Coulomb interaction and the multilevel structure of 
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the QD are considered. However, the spin-flip scattering is only included in the tunneling 
barriers. 

Meanwhile, the significant role of spin-orbit interaction in the QD, which may cause 
the spin rotation of an electron while in the QD, has attracted considerable attentions 
more recently 15 , especially in spin-polarized transport in magnetic nanostructures 16-19 . The 
spin-flip mechanisms in the GaAs-based QD have been investigated in Ref. [15]. Most of the 
theoretical studies 17-19 concentrate on exploring the effect of the intradot spin-flip scattering 
on linear and nonlinear conductances of F-QD-F systems in Kondo regime and a wide variety 
of novel features have been revealed. When the spin-flip scattering strength is of the order 
the Kondo temperature, the original single Kondo peak in the differential conductance is 
split into two peak or three peak structure due to the spin-flip process in the QD 18,19 . Hence, 
it is natural to ask if the intradot spin-flip scattering could induce some novel spectrum of 
tunneling AR conductance in the F-QD-S system. 

In this paper we study the AR tunneling through a F-QD-S hybrid structure by using 
nonequlibrium Green Function method. We mainly emphasize the effect of the spin-flip 
scattering in the QD on linear AR conductance at zero temperature. Until now to our 
acknowledge, there are no theoretical research works to eximine this issue. The spin-flip 
scattering in the QD plays important roles for the AR process of such a F-QD-S system. 
For an isolated QD, it can split one spin-degenerate level of the QD, e d , to two spin-coherent 
levels, e± = Ed ± R, whose states are a superposition of the spin-up and spin-down ones. It 
implicates 17 that incident electrons with up-spin and down-spin from the left F-lead should 
tunnel coherently into the levels split by the intradot spin-flip scattering. In the spin- 
coherent tunneling process, it is expected for the Andreev reflection to bring about some 
novel resonant features of the conductance. We found that the competition between the level 
spliting and the broadening of the split levels arisen from the tunneling coupling, together 
with the spin polarization and the Fermi velocity matching condition, can determine the 
spin-up and spin-down populations of the QD, and further dominates resonant behaviors of 
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the AR conductance of the system. When the spin-flip scattering strength overbears that of 
the tunneling coupling to the leads, the AR conductance versus the gate voltage displays a 
symmetric double peak resonance, and the spin-flip scattering always suppresses the heights 
of the double peaks. However, for a weak spin-flip scattering process in the QD, it only leads 
to a single peak resonance of the AR conductance. In this the spin-flip scattering 

strength increases, the height of the conductance peak may be first increased gradually and 
then dropped, depending the matching condition of the Fermi velocity. 

2. The model and formulas 

Consider resonant AR tunneling through a QD with the intradot spin-flip scattering 
connected to a F-lead and a S-lead, in which only one spin degenerate energy level is included 
and the Coulomb repulsion is ignored for simplicity. The spin quantization axis of the F- 
QD-S system is taken as the direction of the F-lead magnetization, along z-axis. The model 
is shown schematically in Fig.l. The Hamiltonian of the system under consideration, can 
be written as 

H = H F + H s + H dot + H T (1) 

with 

H F = J2( £ k« + ° M )flh« (2) 

k,a 

H S = J2 £ P* s i° s P° + J2( A * s Iv s -pI + As pV s -pi) ( 3 ) 

p,a p 

H dot = ]T e d dld a + + 4d T ) (4) 

<j 

H T = ECJWlA + H.C.) + £(T pCTS t A + H.C.) (5) 

where Hp and Hs are the Hamiltonians for the F-lead and the S-lead, respectively. Under 
mean-field approximation, the F-lead is characterized by a molecular magnetic moment M. 
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The title angle between molecular magnetic moment and the F-QD interface is chosen to 
zero. The BCS Hamiltonian is adopted for the S-lead with an order parameter A standing 
for its energy gap. H dot models the QD with single spin degenerate level e d . The spin-flip 
term in the H dot is caused by spin-orbit interaction in the QD 15 ' 17 and R is the spin-flip 
scattering strength. H T describes the tunneling part between the QD and the F-lead and 
the S-lead with the tunneling matrixes T^ a and T pa . The spin conservation is assumed in 
the tunneling barrier processes, which is distinguished from that in Ref. [14]. 

The current flowing into the central region from the left ferromagnet lead can be evaluated 
from the time evaluation of the total electron number in the left lead 13 ' 20 : 



(6) 



(g r ' a r = 
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Here various kinds of Green functions are expressed in 4 x 4 Nambu representation 20 . The 
Green's functions of the electron of the QD can be exactly solved in the terms of Dyson's 
equation, G r ' a = g r ' a + g r,a XT'" G r ' a , in which XT'" is the self-energy due to both the spin-flip 
interaction in the QD and the spin-dependent tunneling couplings to the left and right leads, 
and g r,a is the Green function without both the tunneling coupling and the intradot spin-flip 
scattering: 

/ 

Lu-e d ±i0+ 

u + e d ±iO + 

u-e d ±iO + 

u + e d ±iO + j 

For the F-QD-S system studied, E r,a can be written as J2 r,a = Er + E/o + E«o - Here the 
off-diagonal term of H dot , the intradot spin-flip scattering, is conveniently considered by the 
self-energy 
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The magnetization of the ferromagnet lead is described by introducing the spin polarization 
factor P. Then Y ^ = Tf (l + P) and Tj^ = Tf (l — P) stand for the spin- up and the 
spin-down tunneling coupling strengths to the F-lead, respectively, resulting in the spin- 
dependent linewidths of the QD level. T f = 2-np 1 jT^T k is the spin-averaged coupling 
strength, Tf = |(r^ + Ff±) denoting the tunneling coupling between the QD and the 
F-lead without the internal magnetization. Within the wide bandwidth approximation, the 
self-energy coupling to the F-lead, Y^f is read as T^f- Here Ff can be written as: 
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with P, the spin polarization in F-lead. The self-energy coupling to the S-lead is: 
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where p r s {uS) is the dimensionless BCS density of states: 
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Vlo 2 - A 2 iVA 2 - u 2 
and F s0 = 2Trp™T*T p is the tunneling coupling strength between the QD and the S-lead. p™ 
in r s0 is the normal density of state while the order parameter A = 0. It is convenient to 
introduce the linewidth function matrices for the S-lead: 
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with pf(co) 



u\e(\u 



A)/y/u 2 — A 2 . After a straightforward calculation, we obtain the 



formula of the tunneling current as follows: 



J = Jn + J a 



(13) 



with 




/ 



du)[fi{uj - eV) - f r (uj)\ [G r d Y s G a d T } ] 



i=l,3 



(14) 



and 



J A = \ I dulffa - eV) - f t {uj + eV)\ jj G^T f G a d F ,) 



(15) 



i=l,3 



where and / r are the Fermi-distribution functions in the left and right leads, respectively. 
Jn is the normal tunneling current which is caused by the single quasiparticle or quasihole 
transport, and J a is the Andreev reflection current. In the linear response regime, the 
normal tunneling conductance and the AR conductance are obtained as follows:: 



Since the normal linear conductance is zero, = 0, at zero temperature, only the Andreev 
reflection process contributes to the linear electronic transport of the system. So the total 
conductance G is equivalent to Ga- 

3. The calculated results and discussion 

We constrain ourselves only to discuss linear AR conductance at zero temperature for 
the F-QD-S, in which the energy level of the QD e d) controlled by the gate voltage V g , is 
restricted in the range of the energy gap of the S-lead ( \ed\ < A) and |e d ± R\ < A. In the 
following calculation, both Fermi energies of the F- and S- leads are set to zero, the energy 
gap of the S- lead, A is taken as energy unit and the spin polarization is chosen as P = 0.3. 




(16) 



and 




(17) 
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First we illustrate the effect of the intradot spin-flip scattering on resonant behaviors of 
the AR conductance versus the energy level of the QD, e d . In Fig. 2, let r s0 = 0.1, we plotted 
the AR conductance as a function of Ed in Fig.2(a) with Tf = 0.02 , Fig. 2(b) Tf = 0.1, 
and Fig. 2(c) Tf = 0.2, for some different spin-flip scattering strengths, R = (solid line), 
0.03 (dashed line), 0.05 (dotted line), 0.07 (dot-dashed line), 0.09 (dot-dot-dashed line), and 
0.15 (short dashed line), respectively. In Fig.2(a), Tf < T s0 , it is clearly seen that for a 
weak spin-flip scattering in the range of R = ~ 0.05, the AR conductance displays a 
single peak resonance at the position of Ed = and its amplitude gradually rises till the 
maximum G m = 4e 2 /h, at R m ~ r s0 /2 = 0.05, with the R increasing. This is a perfect 
AR tunneling process. For some stronger spin-flip scatterings R (0.05 ~ 0.06), however, 
the AR conductance displays also a single peak profile at E d = 0, but the amplitude of the 
resonant peak reduces quickly. As the spin-flip scattering further increases, R > 0.06, the 
original single peak of the conductance develops to a well-resolved double peak resonance. 
The peaks appear near by ±R, respectively. Furthermore, the intradot spin-flip scattering 
always suppresses the heights of the resonant bouble peaks. 

Fig. 2(b) presents some curves of the resonant AR conductance for the symmetric tun- 
neling couplings, Tjo = T s0 = 0.1, and other parameters are the same as those in Fig. 2(a). 
Comparing with Fig. 2(a), a strong enough spin-flip scattering R (> 0.08) brings about a 
double peak resonance of the conductance due to the larger broadening of two split levels 
T = (2rj + T s0 ). It is found that the widths of the resonant double peaks enlarges for the 
enhanced broadening of the minority spin, Ff±. In Fig.2(c), Tjo > r s0 , the amplitude of the 
single peak resonance shows a novel feature: as the spin-flip scattering increases, the peak 
amplitude of the resonance is decreased monotonously. It is worth to notice that in the 
presence of the intradot spin-flip scattering, the single peak of the AR conductance exhibits 
characteristic behaviors essentially depending on a effective overlap of the broadening of the 
two split levels. 

To elucidate the evolution of the resonant AR conductance from single peak to double 
peaks, we calculate the magnitude of the AR conductance at Ed = 0, Go, versus the spin- 
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flip scattering strength R. Defining the ratio of the two tunneling coupling strengths r = 
r s o/T/o, the matching condition of the Fermi velocity, TfyFfi = T 2 s0 reads now as P 2 + r 2 = 
1. Fig. 3(a) shows some curves of the AR conductance Go, for a given r s0 = 0.1 and several 
different r /0 = 0.1 (solid line), 0.1/3 (dashed line),0.1/5 (dotted line), 0.1/7 (dot-dashed 
line), 0.1/9 (dot-dot-dashed line). For the case of r > 1, the magnitude of G increases 
firstly to its maximum 4e 2 /h at R m and then drops fastly as the spin-flip scattering strength 
R increases. It should be mentioned that for the r > 1 where the matching of the Fermi 
velocity can never been satisfied, Go should decrease monotonously with the spin polarization 
P increasing and can not reach to the maximum 4e 2 //i 13 ' 14 . Our calculations indicated that 
there must exist, apart from what considered in Ref.[3] and [14], some other mechanisms 
that result in the perfect AR tunneling, G rising to Ae 2 /h. We believe the intradot spin-flip 
scattering may account for it and leaving somewhat discussion in later. For a enough small 
Tfo, Go becomes a very sharp peak, and its maximum position R m approaches very closely 
to R = r s0 /2. This means that if the spin-flip scattering strength slightly deviates from 
r s0 /2, the AR conductance quickly decreases from 4e 2 /h to 0. 

The typical feature showed in Fig. 3 (a) is understood qualitatively as follows. Spin- up and 
spin-down electrons can escape from the QD through the tunneling coupling to the leads, 
which leads to finite resonant broadening of the two spin-coherent split levels (ed = ±-R) by 
an amount V. Here V = T s0 + + = (r s0 + 2Fj ), the linewidth of the split levels, 
delineates the distribution of the density of states (DOS), in which and are spin- 
dependent tunneling rates to the F-side. r s0 is spin-independent tunneling rate to the S-side. 
When R < R m (— r s0 /2), the linewidths of the two split levels are overlapped effectively 
at Ed = 0, so that the AR conductance versus Ed behaves as single peak resonance. In 
this situation, the AR conductance Go, is enhanced with increasing R, because the intradot 
spin-flip scattering not only shift the level position of the QD from e d = to Ed = ±-R, 
but also change the spin- up and spin-down distribution of the DOS for the split levels 17 . 
Since the minority spin population near the Fermi energy determines behaviors of the AR 
tunneling, the spin-flip scattering turns effectively the majority spin carriers to minority 
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ones near Ed = resulting in Go to rises till its maximum 4e 2 /h, at R m , in which spin- 
up and spin-down carriers from the F-lead completely form pairs into the S-lead. When 
R > T s0 /2 + Tfo > R m , the two split levels have been shifted sufficiently away from each 
other leaving a vanishing spin-dependent DOS at Ed = 0. Therefore G drops quickly to 
zero and it should appear a deep valley in the resonant conductance curve. This implies 
that the AR conductance has developed into a well-resolved double peak resonance shown 
in Fig. 2(a). Fig. 3(b) presents the curves of the AR conductance, G versus R with a fixed 
Tfo = 0.1 and several different r s0 , for r > 1. The peaks exist at a larger R m than that 
in Fig. 3(a) owing to the stronger tunneling coupling rate r s0 , but their patterns are very 
analogous to each other due to the same spin minority Y^. 

In Fig. 4(a), we plotted G as a function of the spin-flip scattering strength R with a 
fixed r s0 = 0.1 and several different Tf = 0.1 (solid line), 0.3 (dashed line), 0.5 (dotted 
line), 0.7 (dot-dashed line), 0.9 (dot-dot-dashed line). This is the situation of r < 1, and 
the magnitude of G decreases monotonously with R increasing. Since the linewidths 
and Tfi are much larger than r s0 , the spin-up and spin-down DOS are compatively low. 
With the increasing of spin-flip scattering, the minority spin occupation reduces gradually 
at Ed = 0. Simultaneously, majority spin carriers can scarcely turn to minority ones near 
Ed = because of the requirement of the energy conversation. As a result, the magnitude 
of G decreases monotonously with the R increasing. In Fig. 4(b), we present some curvers 
of G for the case of r < 1 with a fixed T f = 0.1 , but for several different r s0 - Similar 
features, but a even faster drop of Go with R as in Fig.4(a), have been indicated. As is 
well-kown, r s0 describes the probability that two electrons in the QD tunnel into the S-lead 
and form a Cooper pair. So the weaker the T s0 , the less the probability, and the faster does 
the G decrease to zero as the R increases. 

4. Conclusion 

In summary, we have studied the spin-dependent AR tunneling through a F-QD-S struc- 
ture by using nonequilibrium Green function method. We found that the coherent spin-flip 
scattering in the QD plays important roles in the spin-dependent AR tunneling through the 
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F-QD-S system. The observed single or double peak resonant behaviors of the AR conduc- 
tance, versus the gate voltage, is a consequence of the competition between the spin-flip 
scattering and the resonant broadenings of the two split levels due to the tunneling coupling 
to the leads. When the spin-flip scattering strength in the QD is smaller than the broaden- 
ings of the split levels, the AR conductance exhibites a single peak resonances. In this case, 
as the spin-flip scattering strength increases, the height of the single peak conductance may 
be first increased gradually and then deduced dropped quickly. However, when the spin-flip 
scattering induced spliting of the spin-degenerate level overbears the broadening of the split 
levels, the AR conductance appears as a symmetric double peak resonance, for which a 
novel structure in the tunneling spectrum of the AR conductance is predicted to appear. 
We expect the present results may have practical applications in the field of spintronics. 
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FIGURES 



Fig.l. The quantum dot with intradot spin-orbit interaction is coupled to a ferromagnet 
and a superconductor. A level of the QD is split into two spin coherent levels by the 
spin-flip interaction. 

Fig.2. The resonant curves of the AR conductance versus the energy level of the QD, e d) 
with parameters P = 0.3, r s0 = 0.1 and several spin-flip scattering strengths R = 
(solid line), 0.03 (dashed line), 0.05 (dotted line), 0.07 (dot-dashed line), 0.09 (dot-dot- 
dashed line), and 0.15 (short dashed line) for three different spin-averaged tunneling 
couplings to the F-lead: ( a) Tf = 0.02, Tf < T s0 , (b) Tf = 0.1, Tf = T s0 , and 
Tfo = 0.2, Tf > T s0 . 

Fig. 3. The AR conductance at Ed = 0, Go versus the R with a parameter P = 0.3 (a). 
Tf < r s0 and r s0 = 0.1, the curves of the conductance for some different Tf = 0.1 
(solid line), 0.1/3 (dashed line), 0.1/5 (dotted line), 0.1/7 (dot-dashed line), 0.1/9 
(dot-dot-dashed line), (b). T/o < r s0 and Tf = 0.1, the curves of the conductance 
for T s0 = 0.1 (solid line), 0.3 (dashed line), 0.5 (dotted line), 0.7 (dot-dashed line), 0.9 
(dot-dot-dashed line). 

Fig. 4 The G versus R with a parameter P = 0.3, T/o > r s0 (a). T s0 = 0.1, the curves of 
the conductance for some different Tf = 0.1 (solid line), 0.3 (dashed line), 0.5 (dotted 
line), 0.7 (dot-dashed line), 0.9 (dot-dot-dashed line), (b). Tf = 0.1, the curves of the 
conductance for r s0 = 0.1 (solid line), 0.1/3 (dashed line), 0.1/5 (dotted line), 0.1/7 
(dot-dashed line), 0.1/9 (dot-dot-dashed line). 
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